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Abstract 

Spin Hall effects are a collection of phenomena, resulting from spin-orbit 
coupling, in which an electrical current flowing through a sample can lead to 
spin transport in a perpendicular direction and spin accumulation at lateral 
boundaries. These effects, which do not require an applied magnetic field, can 
originate in a variety of intrinsic and extrinsic spin-orbit coupling mechanisms 
and depend on geometry, dimension, impurity scattering, and carrier density of 
the system — making the analysis of these effects a diverse field of research. In 
this article, we give an overview of the theoretical background of the spin Hall 
effects and summarize some of the most important results. First, we explain 
effective spin-orbit Hamiltonians, how they arise from band structure, and 
how they can be understood from symmetry considerations; including intrinsic 
coupling due to bulk inversion or structure asymmetry or due to strain, and 
extrinsic coupling due to impurities. This leads to different mechanisms of 
spin transport: spin precession, skew scattering, and side jump. Then we 
present the kinetic (Boltzmann) equations, which describe the spin-dependent 
distribution function of charge carriers, and the diffusion equation for spin 
polarization density. Next, we define the notion of spin currents and discuss 
their relation to spin polarization. Finally, we explain the electrically induced 
spin effects; namely, spin polarization and currents in bulk and near boundaries 
(the focus of most current theoretical research efforts), and spin injection, as 
well as effects in mesoscopic systems and in edge states. 

Keywords: spin transport, spin Hall effect, spin accumulation, spin 
current, intrinsic spin-orbit coupling, extrinsic spin-orbit coupling, skew 
scattering, side-jump mechanism, semiconductors, graphene 



1 Introduction 

In the simplest version of a spin Hall effect, an electrical current passes through a 
sample with spin-orbit interaction, and induces a spin polarization near the lateral 
edges, with opposite polarization at opposing edges (D'yakonov and Perel', 1971). 
This effect does not require an external magnetic field or magnetic order in the 
equilibrium state before the current is applied. If conductors are connected to the 
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lateral edges, spin currents can be injected into them. Electrical current in a sample 
can also produce a bulk spin polarization, far from the edges, which is not generally 
classified per se as a spin Hall effect, though it is intimately related and is an 
important ingredient of spin Hall calculations. 

Spin Hall effects have received a great deal of theoretical attention recently, 
in part because the subject includes ingredients of spintronics, electrical generation, 
transport, and control of nonequilibrium spin populations, and also because analysis 
has shown that the problem has remarkable subtlety. Theoretical efforts were also 
fueled by recent experimental observations of these effects. 

In the following we consider semiconductors, where the various mechanisms of 
spin-orbit coupling are well-known and can be roughly classified into two categories. 
First, the extrinsic mechanism is only present in the vicinity of impurities and leads 
to spin-dependent scattering, including Mott skew scattering. Second, the intrin- 
sic spin-orbit coupling remains finite away from impurities and can be understood 
as a (often spatially homogeneous) spin-orbit field inherent in the band structure. 
Furthermore, the spin-orbit couplings are strongly symmetry dependent and are 
therefore different for electron and for hole carriers, and are different for two and 
three dimensional systems. Thus, the microscopic processes involved in generating 
a spin Hall effect can depend critically on such system properties. Initially, there 
was hope that spin transport theory can be formulated in terms of universal spin 
currents that would simplify our understanding; however, it turned out that there 
is no such universality. As there is no unique description of the spin Hall effect, we 
should rather refer to it as a set of spin Hall effects. 

There is already a vast amount of theoretical literature available on spin Hall 
effects and on the related spin currents. It is beyond the scope of this article to 
provide a historical overview or to give an explanation of all the theoretical tech- 
niques used. We rather provide an overview of the various mechanisms, explain 
them using intuitive and qualitative physical pictures, and give a summary of some 
key theoretical descriptions and results. 

In contrast, the number of experiments on spin Hall effects is small and an 
overview is straightforward. In the experiment by Kato et al. (2004b), electrical 
currents in three-dimensional n-GaAs layers (2/im thick) induced a spin Hall effect, 
which was optically detected via Kerr microscopy. Measurements in strained samples 
showed little dependence of the effect on the crystal orientation and it was concluded 
that the extrinsic mechanism proposed by D'yakonov and Perel' (1971) was caus- 
ing the spin Hall effect (Kato et al, 2004b). Indeed, the experimental data can be 
described with reasonable accuracy by the extrinsic mechanism in a model based 
on scattering by screened Coulomb impurities (Engel et al, 2005), as well as one 
based on short-range scatterers (Tse and Das Sarma, 2006a). (See Sec. 14.31 ) Similar 
experiments in ZnSe (Stern et al, 2006) were also in agreement with theory, and a 
spin Hall effect was observed at room temperature. In another experiment, Wun- 
derlich et al. (2005) observed a spin Hall effect in two-dimensional layers of p-GaAs 
by detecting the polarization of recombination radiation at the edges of the sample. 
They ascribed this effect to the intrinsic mechanism, which is consistent with the 
magnitude of the observation (Schliemann and Loss, 2005; Nomura et al, 2005). 
Furthermore, in measurements on a two-dimensional electron system in an AlGaAs 
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quantum well, a spin Hall effect was also observed and ascribed to the extrinsic 
mechanism (Sih et ai, 2005). Finally, Valenzuela and Tinkham (2006) observed a 
reciprocal spin Hall effect in Al, where a spin current induced a transverse voltage 
via the extrinsic mechanism (Hirsch, 1999). 

In addition to theoretical works in the traditional sense, there is also a large 
amount of numerical simulations on concrete realizations of disordered systems, e.g., 
on a finite lattice. It has generally been difficult to make direct comparisons between 
numerical simulations and theoretical predictions, in part because theoretical works 
usually assume that the spin-orbit splittings are much less than the Fermi energy, 
while simulations tend to employ larger spin orbit splittings in order to obtain nu- 
merically significant results. See for example (Ando and Tamura, 1992; Sheng, Sheng 
and Ting, 2005; Nikolic et a/., 2005; Li et ai, 2005). 

In Section [2] below, we review the mechanisms for spin-orbit coupling in the 
semiconductors of interest, and we discuss forms of the effective Hamiltonians that 
describe the carriers in various situations. In Section [3l we see how the various effec- 
tive Hamiltonians can influence spin transport and accumulation. We discuss spin 
precession, produced by the intrinsic spin-orbit coupling, as well as skew scattering 
and the so-called side-jump effect, resulting from the extrinsic spin-orbit coupling. 
We introduce Boltzmann-type kinetic equations which can describe spin-transport 
and accumulation in various situations, and we discuss the simpler spin and charge 
diffusion equations which can typically be used in situations where the spin relax- 
ation rate is much slower than momentum relaxation. 

In Section m we discuss explicitly the spin polarization and spin transport arising 
from an electrical current in a semiconductor with spin-orbit coupling. We introduce 
the notion of a spin current and the spin Hall conductivity, and we discuss results 
that have been obtained for these quantities in various situations. We also discuss a 
relation between the spin Hall conductivity and the so-called anomalous Hall effect 
that can result from spin-orbit coupling in a ferromagnet or in semiconductor with 
a spin polarization induced by an external magnetic field. 

Spin currents are not directly observed in experiments, however. If spin relax- 
ation rates are slow, one may expect that spin-currents with a non-zero divergence 
can lead to observable local spin polarizations, in which relaxation of excess spin bal- 
ances the accumulation of spin that is transported into a region by the spin current. 
Furthermore, boundary effects may be important and non-trivial; in the presence of 
an electric current, spin-polarization may be generated directly at a sample bound- 
ary. These issues are discussed in Subsection 14. 5 1 We also discuss briefly mesoscopic 
systems, where all parts of the sample are close to a boundary. 

A different type of spin Hall effect, associated with edge states, has been predicted 
to occur in certain systems that are insulating in the bulk, where the topology of 
the band structure has been altered due to spin orbit coupling. In Section [SI we 
discuss this concept, along with the possibility that such effects may occur and be 
observable in a number of materials. 
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2 Spin-orbit coupling in semiconductors 



For a non-relativistic electron in vacuum, the Dirac equation can be reduced to 
the Pauh equation, describing a two-component spinor and containing the Zeeman 
term. The Pauh equation also contains relativistic corrections — including the spin- 
orbit coupling 

HsO, vac = Avac (T ■ (k X . (1) 

Here, we used Avac = — fi.^/4mQC^ w —3.7 x 10~^ A^, vacuum electron mass mo, 
velocity of light c, and k — p/h. In a semiconductor, we split the total potential V = 
Vcr + V into the periodic crystal potential Vcr and an aperiodic part V, which contains 
the potential due to impurities, confinement, boundaries, and external electrical field. 
One then tries to eliminate the crystal potential as much as possible and to describe 
the charge carriers in terms of the band structure. The simplest systems of this 
sort can be exemplified by electrons in cubic direct-gap semiconductors. Then, the 
minimum of the energy spectrum is usually near the center of the Brillouin zone, and 
the two-fold Kramers degeneracy is the only degeneracy of the spectrum at fc = 0. 
It follows from symmetry arguments, that for slow electrons in such crystals, and 
for slow carriers (electrons and holes) in high symmetry two-dimensional systems, 
the effective single-particle Hamiltonian is 

H,s = ek + V + Hint + H,^t, (2) 
ifi„t = -ib(k).^, (3) 
TJext = A (T • (k X VF) , (4) 

where k is the crystal wave vector relative to the zone center, and we assumed that 
V is only slowly varying on the scale of the lattice constant. Here cr is the vector of 
Pauli matrices for the pseudo spin-i of the Kramer's doublet at A: = 0; it is customary 
called a spin-i system, b (k) is the intrinsic spin-orbit field, with b (k) — — b (— k) 
due to time reversal symmetry. Thus, for a three-dimensional system, b can only 
be present if the inversion symmetry of the host crystal is broken. In the case 
of a two-dimensional system, it is conventional to talk about its two-dimensional 
bandstructure, and to include the confinement potential in and b (k) (instead of 
including it explicitly in V); in this case b can also result from an asymmetry in the 
confinement. 

In contrast, i/cxt does not require broken inversion symmetry of the pure crystal 
or of the structure. It is important to note that A in Equation ^ can be many orders 
of magnitude larger than the vacuum value Avac; this is due to the large spin-orbit 
interaction when the Bloch electrons move close to the nuclei, with velocities that 
are close to relativistic. Both Hi^t and H^xt may be important for the spin Hall 
effect, as we will discuss in this article. We present specific forms of these effective 
Hamiltonians in Sees. 12.3] and 12.41 
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2.1 Band structure of materials with spin-orbit interaction 

We now consider the electron wavefunctions near the forbidden gap of a semiconduc- 
tor. These wave functions are often described by the Kohn-Luttinger k • p method, 
where one expands the Hamihonian in terms of band edge Bloch functions. Here, 
we present a brief overview of this method; more detailed explanations can be found, 
e.g., in Blount's (1962) review article or in the books by Bir and Pikus (1974) and 
by Winkler (2003). 

For a given k, the solutions of the Schrodinger equation are Bloch functions 
e*'' ''u,y',k (r)- Here, v' is the band index and includes the spin degree of freedom. 
The lattice-periodic part u,y'^k (r) of these Bloch functions can be expanded in the 
functions Ujy. k=o(r) = (r|uy^o)i which provide a complete basis when all bands v 
are taken. For semiconductors with a direct gap at the center of the Brillouin 
zone, which we discuss here, one may consider states in close vicinity of k = 0, 
truncate this expansion, and only take the closest bands into account. Therefore, it 
is sufficient to know the matrix elements of the full Hamiltonian H in the truncated 
basis |</)y) = e^^'^\u„fl), i.e., one considers I/'yi^'(k) = {(j)^ \H\(f),,'). More concretely, 
one evaluates 

Hi,^>{k) = £^,5^1,' + 5,y„' H k • {u^,o\-k\u^' ,o) , (5) 

ZniQ niQ 

where i?^ is the energy offset of the band at k = 0, i.e., [p^/2mo -I- Va + (Avac/ft) c • 
(p X Wcr) ] \u„.q) = E^\u,^^o). Further, tt = p-f (Avac/^) (VVcr x cr), and one usually 
approximates tt w p; then the last term of Equation ((S]) is proportional to the 
matrix element of k • p, giving this method its name. Finally, this finite-dimensional 
Hamiltonian H (k) now describes the band structure in terms of a few parameters — 
band offsets and momentum matrix elements of k = Bloch functions — and is well 
suited for analyzing charge carriers. 

Alternatively, one can use a second method and construct a Hamiltonian by 
allowing all contributions (up to some order in k) that are invariant under the 
symmetry operations of the system — the coupling constants are material-dependent 
parameters. For example, when considering the top of the valence band in Fig.[T]and 
in the absence of inversion asymmetry, magnetic field, and strain, the most general 
form up to quadratic terms in k (i.e., in the effective mass approximation), is the 
4x4 Luttinger Hamiltonian, 



mo 



71 + §72 1 y - 73 (k • J)' + (73 - 72) 



(6) 



which is consistent with the cubic symmetry. Here, J — (J^, Jy, Jz) and Jt are the 
angular momentum matrices for spin | , and ji are the material-dependent Luttinger 
parameters. Hi^ describes p-doped Si and Ge; for GaAs, due to broken inversion 
symmetry, terms linear in k arise as well. 

By contrast, using the first method (k • p method) instead, the Hamiltonian 
77yy'(k) [Equation ^] is evaluated directly for the eight bands v (including spin) 
shown in Fig. [U and one arrives at the simplest version of the 8x8 Kane model. It 
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Figure 1: Schematic band structure of a cubic direct gap semiconductor. When 
spin-orbit interaction is disregarded, one finds an s-hke conduction band and a p- 
hke three-fold degenerate valence band. The spin-orbit interaction due to crystal 
potential 14r [entering as V in Equation H])] partially lifts this degeneracy and leads 
to a substantial splitting between the valence bands with total angular momentum 
J = 3/2 (heavy and light holes) and those with J = 1/2 (split-off holes). 

only includes three parameters, namely, the energy gap Eq, the energy of the split- 
off holes Aq, and the matrix element P of the momentum (multiplied by h/mo) 
between s- and p-type states. P is nearly universal for III-V compounds, while 
the other parameters depend on the material. This eight-dimensional description is 
accurate for narrow band materials; for wider gaps it still provides understanding 
at a qualitative level. 

Furthermore, one can then derive an effective, lower-dimensional Hamiltonian by 
block-diagonalizing H^i^i , this is an efficient way to calculate the band structure in 
the vicinity of k = 0. One can do this either exactly or by using time-independent 
degenerate perturbation theory (see Sec. 12. 2p . Considering a particular block, this 
allows estimating the magnitude of the symmetry-allowed terms. Terms that were 
not present before block-diagonalization are called contributions from remote bands. 
For example, using the 8x8 Kane model, one can calculate the parameters 7^ that 
enter the 4x4 Luttinger Hamiltonian [Equation ^] for the top of the valence band; 
because the model is isotropic, one gets 72 =73- To estimate corrections due to the 
cubic symmetry, one needs to take more bands into account. 

In addition to the crystal field, we also wish to consider electric fields that are 
applied externally or that result from charged impurities. Assuming that the cor- 
responding potential V (r) varies slowly on the scale of a lattice constant, we can 
apply the envelope function approximation (EFA), i.e., we replace the plain waves 
gik r slowly varying envelope functions (r). Evaluating matrix elements of 
V in the basis {rlip,!^) — (r)it,y, k=o (r), the main contribution is diagonal with 
respect to the band index i'. However, because k and V (r) do not commute, off- 
diagonal elements (in v) can arise in the Hamiltonian when expanding in k, see e.g. 
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Equation Magnetic fields can be included in a similar way, using the Perierls 
substitution hk = —ihV — (e/c) A. Here, we use e < for electrons and e > for 
holes. Then, the components of k no longer commute, k x k = i (e/Hc) B; this leads 
to Zeeman coupling, which is described by an effective g factor. 

2.2 Effective Hamiltonian 

The k • p method leads to high-dimensional Hamiltonians [Equation ([5])], for exam- 
ple, an 8 X 8 matrix for the Kane model, thus further simplifications are desirable. 
For this, one can use time-independent degenerate perturbation theory and describe 
a subset of states (say, the lowest conduction band states) with an effective Hamilto- 
nian. Lowdin partitioning is a straightforward and convenient method to implement 
such a perturbative expansion (Winkler, 2003); it is also known as Foldy-Wouthuysen 
transformation in the context of the Dirac equation and as Schrieffer-Wolf transfor- 
mation in the context of the Anderson model. The idea is to find a unitary trans- 
formation (i-e., S is anti-Hermitian) such that the transformed Hamiltonian 
e^^He^ is block diagonal, i.e., contains no off-diagonal elements between the states 
we are interested in and any other states. This procedure assumes that the states 
of interest (e.g., the conduction band) are separated from the other states (all other 
bands) by an energy much larger than the Fermi energy. Then, because these off- 
diagonal elements are small, one can eliminate the off-diagonal blocks of e^^ He^ 
order by order (or even exactly). In our example, the transformed Hamiltonian con- 
sists of one 2x2 and one 6x6 block. The smaller block describes the conduction 
band electrons — we can understand the two dimensions as (pseudo-) spin i. At the 
level of wave functions, the periodic part of the electron wave function at a given 
k 7^ is mainly described by the conduction band Bloch function at A: = 0, but also 
contains a small admixture from the valence band Bloch functions. 

2.3 Intrinsic spin-orbit coupling 

One generally distinguishes between intrinsic and extrinsic mechanisms of spin-orbit 
coupling; however, this classification is not unique across the literature. In this 
article, we classify it according to the individual terms of the effective Hamiltonians. 
Namely, we refer to the spin-orbit contributions to the Hamiltonian that depend on 
impurity potentials as extrinsic. The other spin-orbit contributions arise even in 
the absence of impurities and we call them intrinsic — we also call effects resulting 
from these contributions intrinsic, even if we must allow for a small concentration 
of impurities to make the theory of dc transport properties consistent 

For a (pseudo-) spin-i system, the spin-orbit part of the intrinsic one-particle 
Hamiltonian has the general form Hint [Equation ([3])]. In the following, we discuss 
the origin and the functional form of such spin-orbit fields. We focus on such spin-^ 

^This is reminiscent of the definition of an "intrinsic semiconductor," which is so pure that (at 
a sufficiently high temperature) the impurity contribution to the carrier density is negligible. The 
conductivity of such a sample is known as intrinsic conductivity (again, a finite transport lifetime 
T is required to make the conductivity well-defined). At lower temperatures, the carrier density 
mainly results from the impurities and now one refers to extrinsic properties. 
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descriptions, because they are relevant for low-dimensional systems and are the basis 
of most theoretical works. 

We first consider a n-doped bulk (3D) semiconductor and the effective Hamil- 
tonian of conduction band electrons. III-V and II- VI semiconductors lack inversion 
symmetry and are available in two modifications: in cubic zinc blende or in hexag- 
onal wurtzite structure. In zinc blende modification, the bulk inversion asymmetry 
(BIA) leads to the Dresselhaus term 

Hu,3d^Bk,{kl^kl)a^+c.p., (7) 

where ki are along the principal crystal axes. Here, c.p. stands for cyclic permutation 
of all indices, and the symmetrized product of the components ki must be used if 
a magnetic field is applied. The Dresselhaus term originates from bands further 
away than the basic eight bands, and one finds the coupling constant in terms of 
the band parameters; when using the extended 14 x 14 Kane model its numerical 
value is S « 27eVA3, for both GaAs and InAs (Winkler, 2003). However, tight- 
binding calculations and interpretation of weak-localization experiments indicate 
lower values, at least for GaAs (Knap et al, 1996; Krich and Halperin, 2007). 

When the electrons are confined to two dimensions, the expectation value of the 
Dresselhaus term along the confinement direction (that we always assume to be along 
[001]) should be taken, (-ffD,3d)- While (fc^) — 0, we see that the terms in (fc^) « 
(Tr/d)^ are large for small confinement width d, thus the main BIA contribution 
becomes 

Hp = (3 (k^cr^ - kyay) , (8) 

with (3 « (tt/c?)^. In addition to the fc-linear term in Equation ([5]), there is also 
a fc'^-term, 

Hr,^2d = Bk^ky {kyax - k^ay) , (9) 

which is small compared to Hp in the strong confinement (low carrier density) limit 
ir/d ^ fcp, where fcp is the Fermi wave vector. Additionally, a spin-orbit coupling 
term arises if the confinement potential V (z) along the z-direction is not symmetric, 
i.e., if there is a structure inversion asymmetry (SIA). Equation (|4]) provides an 
explicit connection between such a potential and spin-orbit coupling. Taking the 
expectation value (i?cxt) along the z-direction and noting that the only contribution 
of the confinement field is oc {VzV), one finds the Rashba Hamiltonian, 

Ha ^ a {kyax ~ k^cFy) , (10) 

corresponding to b (k) = 2a z x k. More generally, for spinors with — ±1/2, 
Ha is the only fc-linear invariant of the group Coov that takes into account the 
confinement potential V{z) but disregards the discrete symmetry of the crystal. 
The magnitude of the coupling constant a depends on the confining potential and it 
can be modified by applying an additional field via external gates. It also defines the 
spin-precession wave vector ka = am/h^ . Finally, such a term Ha is also present for 
three-dimensional electrons in systems of hexagonal wurtzite structure (or in cubic 
systems with strain, see Sec. 12. 5|) . 
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Next, we consider a p-doped three dimensional semiconductor, i.e., the J = 3/2 
valence band, described by the four-dimensional Luttinger Hamiltonian. Remote 
bands lead to a BIA contribution to the Hamiltonian, which is given by Equation 
([7]) after replacing ai by the angular momentum matrices Ji for spin | and using 
a different coupling constant. If the system is reduced to two dimensions, size 
quantization lifts the fourfold degeneracy at k = and creates heavy-hole (HH) 
bands with — ±3/2, and light-hole (LH) bands with = ±1/2 (for confinement 
along the [001] axis). Usually the HH bands are higher in energy, thus for small 
doping it is sufficient to consider only them. For spinors with Jz = ±3/2, the 
only invariant of the group Coov (again, we do not discuss invariants of discrete 
symmetries here) respecting time reversal symmetry is 

Ha,h = iah {kl<T+ - k\cj^) , (11) 

where a± = ± iay for any a. As distinct from Ha [Equation (jlOp ], the Rashba 
Hamiltonian for heavy holes is cubic in fc, as it was discussed in (Winkler, 2000; 
Schliemann and Loss, 2005). 



2.4 Extrinsic spin-orbit coupling 

Electric fields due to impurities lead to extrinsic contributions to the spin-orbit 
coupling. Externally applied electrical fields lead to analogous contributions. To 
derive the dominant extrinsic term, it is sufficient to restrict ourselves to the simplest 
8x8 Kane Hamiltonian; higher bands will give rise to small corrections. Using third- 
order perturbation theory and for conduction band electrons, we find -ffcxt as given 
in Equation (j4]), with (Nozieres and Lewiner, 1973; Winkler, 2003) 

^0 (Eo + Ao) _ 

and where V is the potential due to impurities and an externally applied field. It is 
noteworthy that Equation 111 has the same analytical form as the vacuum spin-orbit 
coupling [Equation ([1])]; this is because both the Dirac equation and the simplest 
Kane Hamiltonian have spherical symmetry and because both the Pauli equation 
and Equation ([4]) are obtained in a low-energy expansion. However, for Ao > the 
couping constant A has the opposite sign as in vacuum. 

One finds A w 5.3 for GaAs and A w 120 for InAs, i.e., spin-orbit coupling 
in n-GaAs is by six orders of magnitude stronger than in vacuum, and even larger for 
InAs due to its smaller gap. This enhancement of spin-orbit coupling is critical for 
developing large extrinsic spin currents. Furthermore, for a two-dimensional system, 
when considering V as averaged along the z-direction, both \7V and k are in-plane, 
thus we have 77cxt,e = A Cz (k x VF)^. 

For a 3D hole system, we consider the J — 3/2 valence band. Then, the dominant 
extrinsic spin-orbit term in third order perturbation theory describing the valence 
band states is 

-ffext,v = AvJ-(kx Vy), (13) 
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with Av = -P'^/Wl, i.e., for GaAs K ~ -15 (Winkler, 2003). and has to be 
added to the Luttinger Hamiltonian Hi, [Equation ([6])]. When considering a two- 
dimensional hole system with HH-LH splitting, we can restrict Equation (|13p to the 
heavy holes states, where Jz — ±3/2. Expressing this two-dimensional subspace in 
terms of a pseudo-spin ^ leads to 

i?cxt,v = -:^ (kx Vy),a.. (14) 

Thus, the extrinsic spin-orbit interaction for two-dimensional heavy hole states has 
the same form as for two-dimensional electrons. 

Finally we point out that extrinsic spin-orbit coupling arises because the the 
long range Coulomb potential of the impurities does not commute with the intrinsic 
Hamiltonian of the hosting crystal. The extrinsic Hamiltonian i?cxt [Equations (jl]) 
and (fT3|) ] is obtained in the framework of the EFA, which disregards short-range con- 
tributions to the spin-orbit coupling arising from the chemical properties of dopants. 
This is why the coupling A depends only on the parameters of the perfect crystal 
lattice. 



2.5 Strain 

Non-hydrostatic strain reduces the symmetry of the system and in this way leads to 
additional spin-orbit terms in the Hamiltonian. In third order perturbation theory 
of the Kane Hamiltonian, the effective conduction-band Hamiltonian due to strain 
is dominated by 

= Si^'o'S + Ac) " • " ''^ - ^^3- • (k X s.), (15) 

where £s = (Syz, £xz, £xy) describes the shear strain. Here, C2 is the interband 
deformation-potential constant that arises in noncentrosymmetric semiconductors 
(Winkler, 2003; Trebin et al, 1979). Note that Pikus and Titkov (1984) as weU as 
Ivchenko and Pikus (1997) use the opposite sign in the definition of this constant, 

PT /IP 

C2 = — C2. Further, if a shear is applied such that only Exy 7^ 0, Equation p5|) 
has the same form as the Rashba Hamiltonian [Eq. (fTU]) ]. 

For three-dimensional J — 3/2 valence band states, the main strain contribution 
is (Pikus and Titkov, 1984) 

iy,^v = -^j.(kx£,). (16) 

Note that when the system is confined to two dimensions. Equation (I16p 
inplies a /c-linear spin-orbit contribution for heavy hole states due to strain, 
{C2P/E0) (k X es)^<Jz- This linear term arises due to the low symmetry of the 
strained material, in contrast to the Rashba term [cf. Equation (fTTj) ]. which is 
dominated by terms cubic in k and where the fc-linear terms are numerically small 
(Winkler, 2003). 
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2.6 Anomalous velocity and coordinate 



As a consequence of the spin-orbit interaction, velocity and coordinate operators are 
modified and become spin dependent — this will be important when considering cur- 
rents. When an effective Hamiltonian is derived in perturbation theory, as explained 
in Sec. 12.21 a unitary transformation e^'^ is applied. Thus, the coordinate operator 
r = i (d/dk.) is also transformed, r i-^ f = e^'^re'^ = r+5r and we call Sr the anoma- 
lous coordinate. In particular, because S couples to spin, Sr is spin-dependent. This 
correction 6r is known as the Yafet term (Yafet, 1963); also, it can be expressed 
in terms of a Berry connection, (u;/', k |«Vk| k) (Kohmoto, 1985; Sundaram and 
Niu, 1999; Nagaosa, 2006). In perturbation theory, one finds 



for conduction band electrons and for J = 3/2 heavy hole states, respectively. Note 
that coordinate operators no longer commute, r x r = iXcr and r x r = zAvJ, 
resp. Finally, <5rso leads to an extra term in the equations of motion that can be 
understood as anomalous velocity (Blount, 1962). 

Formally, we can derive the anomalous velocity similarly to the coordinate, 
namely v i— > e~"^ve'^ = vq -I- Sv, where Sv is the anomalous velocity and, for a 
parabolic band, vq = Kk/m*. Alternatively, one can obtain the velocity operator 
from the Heisenberg equation, v = (i/h) [H, f]. For H = Hq + Hso, where Hgo con- 
tains the (small) spin-orbit couphng, we get v = vo+{i/h) [Hso, r]-t-(i/ft) [Hq, Srso]- 
Thus, Hso leads to an anomalous velocity because it does not commute with the 
unperturbed coordinate r, and, additionally, the contribution from the anomalous 
coordinate Srgo should also be taken into account, as it can be significant. 

When the impurity potential V is included, the above argument remains the 
same, but now Hso contains the extrinsic contribution as well. Note that for ex- 
trinsic spin-orbit Hso = ^cxt,c [Equation (jj])] the commutator [-ffext.e, r] and the 
anomalous coordinate Srso give equal contributions to Sv (Nozieres and Lewiner, 



3 Mechanisms of spin transport 

We now address how the microscopic mechnism of spin-orbit interaction, given by 
effective Hamiltonians, influences spin transport and accumulation. In the following, 
we assume a non-interacting system in the absence of a magnetic field. Because 
we ignore electron-electron interaction, we do not consider the spin-drag effect here 
(Hankiewicz and Vignale, 2006), which can lead to a suppression of spin transport at 
high temperatures, and suppression of spin-relaxation (Glazov and Ivchenko, 2002). 
We also restrict ourselves to Boltzmann transport and do not discuss the hopping 
regime (Entin-Wohlman et ai, 2005). 



<5rso,c = A(cr X k) , 
^rso.v = Av (J X k) , 



(17) 
(18) 



1973). 
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3.1 Intrinsic: spin precession 

In a system with weak intrinsic spin-orbit coupling iJjnt [Equation ([3])], consider a 
carrier with spin ahgned along the spin-orbit field b (k) . When an electrical field 
E = Ex is applied, the particle is accelerated: k = eE/h in lowest order in spin 
orbit-interaction; and its spin-orbit field changes: b = {dh / dk^) eE / h. For a small 
acceleration, the spin follows adiabatically the direction of b (k) . Additionaly, there 
will be a non-adiabatic correction that can be derived as follows. Say, the direction 
of b rotates in the xy plane (as it is the case for Rashba interaction Ha)- Because 
the rotation frequency oj is the component of b/6 perpendicular to b, it is w = 
{h xh)z/h'^ ■ In the co-rotating frame, there is a field h along the x-axis and a field 
IvjJ along the z-axis. As we are interested in the next-to-adiabatic correction, we 
assume that w changes slowly and that the spin remains aligned along the total 
field in the rotating frame, thus it has a component ~ (^/2) hui/b. Therefore the 
first non-adiabatic correction to the spin is Ss (k) = h^{h x b)^/2&'^. In particular, 
the electrical field drives this spin precession (via b), leading to such out-of-plane 
component Ss (k), which could be important in spin transport (Sinova et al, 2004), 
see Sec. 14.31 However, when considering dc properties, one must be careful and also 
allow for impurities that decelerate the carriers to reach a steady state. In particular, 
if b (k) is linear in k, it turns out that the deceleration at impurities cancels this 
spin precession, see Sec. 14.31 

3.2 Extrinsic: skew scattering 

When a carrier scatters at an impurity potential V, because of the extrinsic spin- 
orbit interaction [Equations ([4]) and ([14])] the scattering cross section depends on 
the spin state (Smit, 1958), see Fig. [2l This effect is known as Mott skew scattering 
(Mott and Massey, 1965) and was originally considered for high-energy electrons 
that are elastically scattered by an atom and that are described by the vacuum 
Hamiltonian, Equation ([T]). Skew scattering does not appear in the first order Born 
approximation, thus is at least of the order V^. For band electrons, skew scattering 
was originally considered as the origin of the anomalous Hall effect, see Sec. 14.41 As 
applied to spin Hall effct, the relevance of this extrinsic mechanism was recognized 
early on (D'yakonov and Perel', 1971; Hirsch, 1999; Zhang, 2000). 

3.3 Extrinsic: side jump mechanism 

The side jump mechanism (Berger, 1970) describes the lateral displacement of the 
wave function during the scattering event. (Such a displacement does not modify 
the skew scattering cross section introduced in Sec. 13. 2[ because it does not change 
the scattering angle measured at large distances.) The side-jump contribution is 
obtained when the anomalous velocity 5v (see Sec. 12. 6|) is integrated over the du- 
ration of the scattering process. As indicated in Sec. 12.61 the anomalous coordinate 
for electrons [Equation ^T7} ] leads to an anomalous velocity contribution A(cr x k); 
and there is an equal term due to {i/fi)[Hcxt,c, r]. For impurity scattering with 
momentum transfer Sk., this results in a total lateral displacement 2A(<t x (5k) [and 
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Figure 2: Spin-dependent scattering of electrons at an attractive impurity. We 
show the classical trajectories (solid lines), for a screened Coulomb potential and for 
strongly exaggerated extrinsic spin-orbit coupling [using Equation ([4]) with A > 0] 
and with spin quantization axis perpendicular to the plane. The skew-scattering 
current results from different scattering angles for spin-]" and spin- J, electrons and 
leads to a positive spin Hall conductivity, CTgg = — jgg y/E^ > 0. The dashed lines 
show the horizontal displacement due to the side jump effect, contributing to the 
spin current with opposite sign. 



anologousely for holes, using Equation (jlSpj. When also the effect of the anomalous 
velocity due to an applied electrical field is considered, the side-jump contribution 
to spin-transport becomes a subtle issue, for a detailed analysis and intuitive de- 
scription see (Nozieres and Lewiner, 1973). Because the side jump mechanism is not 
contained in the Boltzmann approach, in such a framework it needs to be evaluated 
separately. The side jump contribution can be found using the Kubo formula and 
diagrammatic approaches; see (Lewiner et al, 1973; Crepieux and Bruno, 2001; Tse 
and Das Sarma, 2006a). 



3.4 Kinetic equation 

We consider the Hamiltonian Hcs [Equation ([2])] containing both intrinsic (i/int) and 
extrinsic (.Hext) spin-orbit contributions, and with V that describes the electrical 
field E and the impurity potential VJ. For a homogeneous system of non- interacting 
particles, one can derive the kinetic equation (Khaetskii, 2006; Shytov et al, 2006) 

\ / coll. 

i.e., a spin-dependent Boltzmann equation, where the distribution function is written 
as a 2 X 2 spin matrix / = /o(k) -t- ^?^(k) 1 -I- f (k) • cr, with equilibrium distribution 
function /q. Here, n is the excess particle density and f describes the spin polar- 
ization density. Formally, the Boltzmann equation is obtained by an expansion in 
l/fcp^, where ^ is the mean free path, i.e., it is applicable for dilute impurities. Tra- 
ditionally, the Boltzman equation describes the distribution function n(k), which is 
the probability density of a state k to be occupied — in contrast, here it describes 
/(k), which corresponds to the 2x2 density matrix for a spin-i particle. 
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All terms on the l.h.s. of Equation arise in the absence of impurities. The 
first term is the derivative of / with respect to its explicit time-dependence. The 
second term describes the spin precession; it is obtained from the Heisenberg equa- 
tion, f • <T = f • [ — ib (k) • (T, cr] . The third term is the driving term due to 
the electrical field, given in lowest order of E. Finally, for inhomogeneous particle 
and spin distributions, the term v • V/ has to be added to the l.h.s. 

The r.h.s. of the Boltzmann equation [Equation (fT9|) ] is the collision term, sym- 
boUically given by 



(dm] 

\ dt / ^ dn 

\ J coU k': e'=e 



E ^ [/(k')-/(k)l, (20) 



coU 

Here, rii is the impurity density and we only consider elastic scattering k — s- k' and 
k' k. The scattering cross section tensor da/dfl and the summation over final 
states k' are spin-dependent because (i) the extrinsic interaction Hext leads to skew- 
scattering, (ii) the intrinsic spin-orbit Hamiltonian Hmt induces a spin-dependent 
density of states (DOS), and (iii) Hint causes spin-dependent momentum transfer — 
in general, the spin dependence of scattering is rather complex. The description 
of scattering is simplified for weak spin orbit coupling, as the collision term can 
be expanded in spin-orbit coupling, and we can discuss the individual corrections 
seperately. Note that the Boltzmann equation does not include the side-jump effect 
(cf. Sec. 13. 3p . which is of higher order in l/kpi. 

When considering only spin-orbit coupling due to i?cxt , the collision term includ- 
ing skew scattering for a central symmetric impurity potential is (Engel et al. , 2005) 



coU 



n,v, I dn{k')\l{d) /(k')-/(k) +i/(z9)5(^?)cr.m[n(k)+n(k') 



(^) 

(21) 

where m = k' x k/ |k' x k| is the unit vector normal to the scattering plane and 
"ff = ^^kk' is the angle between k' and k. The coefficient !{'&) is the spin-independent 
part of the scattering cross section, while S{'d) is the so-called Sherman function 
(Mott and Massey, 1965; Motz et al, 1964; Huang et al, 2003), describing the 
polarization of outgoing particles (which is normal to the scattering plane) scattered 
into direction k from an unpolarized incoming beam of momentum k'. Note that 
/ S, as mentioned earlier in Sec. 13.21 vanishes in first-order Born approximation; the 
lowest term is cx . Also, S is proportional to spin-orbit coupling, this is why in the 
second term in Equation (j2ip . only the spin-independent part of /, in, is retained. 

So far, we considered the distribution function / (k) as density in k-space. In 
the presence of intrinsic spin-orbit interaction, the energy spectrum contains two 
branches: for a given k, there are two energies, split by the intrinsic field b (k). Thus, 
for elastic scattering, energy is conserved but |k| is not. It is now more convenient 
to choose a distribution as function of energy e and direction ip in k-space. Namely, 
the distributions $c (</5, e) and $ (</?, e) are derived from the distributions n (k) and 
f (k), and can again be writen as matrix ^{ip, e). Note that contains the spin- 
dependent DOS. We now consider a two-dimensional system, k — (k cos ip, ksinip), 
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and assume that the spectrum Ck in the absence of spin-orbit interaction is isotropic 
[cf. Equation ([2])], and we define such that e^^ = e and define — f-'kj^- For 
E = E'k and for b ^ E-p, the kinetic equation becomes (Shytov et at, 2006) 



'dt' 



b X $ 



b X 



dh 

9fc 



eE 



2 de 



kx Id 
— H 

h 2h?v^ dip 



(b • cr sin (p) 




where fo is the Fermi distribution function and b is evaluated for |k| = k^. 

Now, considering only iJint, the collision integral can be found in a Golden Rule 
approximation (Shytov et at, 2006), 



d^iv, e) 
dt 



coll 



dip' K{^) ^{ip') - $(^) 



+ / V rr • [M (ip, (^') - M (^', (^) <i>e (v?)] . (23) 

Jo 

Here, the first term describes the spin-independent scattering, with 
K{d) = K{(p'-ip) = W {q)kj2TTh^v, and q = 2fc, sin (|79| /2). The factor 
^{q) = ( |M(q)|^ ) does not depend on the direction of the momentum transfer q 
because the problem is isotropic (while of course the individual scattering event is 
anisotropic, i.e., depends on the scattering angle i?). This spin-independent term 
coincides with first term of Equation (|2ip . with K {■&) — fiiVel {"&), and should only 
be included once. The second term in Equation (|23p is given in first order in the 
intrinsic spin-orbit interaction b and contains the kernel (Shytov et al, 2006) 



keh (ip) 



b(^)+b(^0 dK j,) , 

AhKv, dd 2' ^ ' 



where the first term results from the spin-dependent DOS of the outgoing wave. 
The second term in M arises, because for a given energy e, |k| depends on the spin 
state. Thus the incoming and outgoing states can have different momenta, leading 
to spin-dependent corrections to q. 

For a very smooth scattering potential such that typically q < h/vp, the spin 
motion is adiabatic and should be treated differently (Govorov et al, 2004; Khaetskii, 
2006). 



coll. 



3.5 Diffusion equation 

A spin diffusion equation can be derived, starting from the Boltzmann equation, for 
a dirty system when the spin relaxation time Tg is much longer than the momentum 
relaxation time r, i.e., Tg ^ r. It describes the carrier density (r) and spin polar- 
ization density s (r); say, is the excess density of particles polarized along z. For 
conduction band electrons, the (pseudo-) spin density is S = {h/2)s. The diffusion 
equation is simpler to solve than the kinetic equation (jl9p . as the dependence on k is 
integrated out. Also, it is usually sufficient to know s, because it an experimentally 
accessible quantity, while /(k) is not directly accessible, cf. Sees. 14.2 1 and 14.51 below. 
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For a two-dimensional system with Rashba spin-orbit interaction Ha [Equation PU]) ]. 
the diffusion equation is (Burkov et al, 2004; Mishchenko et ai, 2004) 

N^DV^{N + poVE)+Tsciyxs)^, (25) 
- D\/\s^ - rrh, + T^s [(z x V) x s], + r^c (z x V), {N + pqVe) , (26) 

with diffusion constant D = jfpT, anisotropic Dyakonov-Perel (1972) spin relax- 
ation rates t^^ = Ty^ = rj^ = 2T{akF/h)'^ and t^^ ~ 2rJ^, spin-charge cou- 
phng Fsc = —2a (afcpr)^ /h^, spin-spin coupling Fss = AaEFr/fi^, density of states 
po = m/nh?, and potential energy Ve of a carrier in the electrical field. The charge 
current is 

r = -DV{N + poVE)+T,,ixs. (27) 

Further, Mal'shukov et al. (2005) derived diffusion equations for two-dimensional 
electrons with the Dresselhaus Hamiltonian i?D,2d [Equation Q]. 

The boundary conditions of the diffusion equation for a system with spin-orbit 
interaction are not trivial and one expects that they depend on the microscopic prop- 
erties of the boundaries. A number of papers have been written about the bound- 
ary conditions corresponding to various physical circumstances and have partly 
clarified this issue (Govorov et al, 2004; Mal'shukov et al, 2005; Adagideli and 
Bauer, 2005; Galitski et al, 2006; Bleibaum, 2006; Tserkovnyak et al, 2006). De- 
pending on the particular boundary condition, there may or may not be an Sz spin 
accumulation near the boundary of a 2D system, see Sec. 14.51 Somewhat related 
to these question, Shekhter et al. (2005) considered the boundary between a diffu- 
sive and ballistic system and allowed for spin-dependent scattering at the boundary, 
resulting from a spatially dependent Rashba Hamiltonian Ha. 



4 Electrically induced spin polarization and spin 
transport 

4.1 Spin current and spin Hall conductivity 

In this article, we define the spin current in a homogeneous system with density n 
as ^ 

where (•) is the expectation value of single-particle operators and with (1) = 1. Thus, 
the spin current is defined as the difference of the particle currents densities (mea- 
sured in numbers of particles) for carriers with opposite spins. This is in accordance 
with many studies (Murakami et al, 2003; Sinova et al, 2004; Sinova et al, 2006), 
where a definition as in Equation (^5)) was chosen, up to numerical prefactors. In 
many definitions, an additional prefactor of ^ is used, which results from h/2 angular 
momentum per electron spin and setting h — 1. With the same argument, for the 
HH band, sometimes a prefactor 3/2 is used; but sometimes only a factor i is used 
to have the same definition of j'' for electrons and holes. Furthermore, the r.h.s. of 
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Equation (pS)) is sometimes multiplied by the charge e to obtain the same units for 
charge and for spin currents. In particular, this means that the sign of the definition 
of j may change if e < for electrons is taken. 
Next, we define the spin Hall conductivity 

= ^31/ E,, (29) 

where jy is the spin current density resulting from a small applied electrical field ■ 
The negative sign in Equation (j29p results from writing a formal similar definition 
for cr^^ as for the charge conductivity (Jxy', however, sometimes a definition with an 
opposite sign for is used. 

These various prefactors are only some technicality — the main question is 
whether defining as in Equation (|28p makes sense. To describe spin transport it 
sounds attractive to find a scheme similar to the charge transport theory. Because 
of charge conservation, charge densities p"^ and charge currents j'^ satisfy the 
continuity equation p*^ + divj'^ = 0. For spin transport, we can consider the spin 
density Si instead of p'^. Mott's (1936) two-channel model of electron transport in 
ferromagnetic metals is based on independent and conserved currents of up- and 
down-spin electrons, and S and jl. obey a continuity equation. The Definition (j28p 
is the natural generalization of Mott's model; however, spin-orbit coupling violates 
spin conservation, and the continuity equation for spin densities and currents does 
not hold. In this article, we will still use Equation (|28p as definition of the spin 
current, as it is widely used, but we remain aware of its limitations. Despite the 
fact that it cannot be directly related to spin accumulation, it is a useful model 
quantity to compare the effect of different spin-orbit coupling mechanisms. While 
the continuity equation does not hold, one can, for a concrete Hamiltonian, evaluate 
source terms arising on the r.h.s. (Burkov et al, 2004; Erlingsson et al, 2005), 
which is often termed as torque (Culcer et al, 2004). 

Other definitions of spin currents have also been proposed. Zhang and Yang 
(2005) analyzed the current of the total angular momentum + Sz and argued 
that it vanishes for the Rashba Hamiltonian Ha in the absence of impurities (due 
to the rotational invariance of Ha) and that thus the impurity scattering would 
determine angular momentum currents. Shi et al. (2006) discussed spin currents, 
introducing a definition that is not proportional to our j^, but is given as time- 
derivative of the "spin displacement" Si (r) rk ■ A somewhat related definition was 
used by Bryksin and Kleinert (2006), who found that such spin currents diverge 
when the frequency a; — s- 0. 

4.2 Spin polarization 

In experiments, the spin polarization can be detected optically (Meier and Za- 
kharchcnya, 1984). Electrically induced polarizations were inferred from measure- 
ments of the Kerr rotation, where the polarization of a linearly polarized beam of 
light rotates when the beam is reflected at a spin-polarized sample (Kato et al, 
2004b; Sih et al, 2005). Alternatively, the circular polarization of the recombina- 
tion light at a p-n junction can be used to determine the initial polarization of the 
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carriers (Wunderlich et al, 2005). Finally, the inverse effect, the photo- galvanic ef- 
fect, can be observed, where a spin polarization is produced by polarized light and 
the induced electrical current is detected (Ganichev and Prettl, 2003). 

In the bulk of a two- or three-dimensional sample, spin polariziations arise be- 
cause an electrical field shifts the Fermi sea; (k) = eEr/h for small E and with trans- 
port lifetime r. This implies that due to intrinsic terms, there is on average a finite 
spin-orbit field, (b (k)). This leads to a bulk spin polarization, which, in simple cases, 
is aligned is along (b) (Ivchenko and Pikus, 1978; Vas'ko and Prima, 1979; Levitov 
et al, 1985; Edelstein, 1990; Aronov et al, 1991). Such a polarization was observed 
experimentally in two-dimensional GaAs hole systems: Silov et al. (2004) used a 
(OOl)-surface and detected the polarization of the photoluminescense from the side 
of the cleaved sample, while Ganichev et al. (2006) used samples of low crystallo- 
graphic symmetry and detected the polarization along the growth direction. Fur- 
thermore, in the presence of anisotropic scattering (see Sec. 13. 4p . a magnetic field B 
can lead to polarization perpendicular to both (b) and B (Engel et al. , 2007) — such 
a perpendicular polarization was already observed by Kato et al. (2004a). 

4.3 Spin currents in bulk 

The bulk spin current was analyzed for many different spin-orbit Hamiltonians. 
For a two-dimensional electron system, the intrinsic effect of the Rashba coupling 
Ha lead to some debates. Because jl. is invariant under time reversal, it is allowed 
to be finite in equilibrum. Indeed, such equilibrum spin currents are predicted 
theoretically, however, they are of order and usually small (Rashba, 2003). 

Now we discuss bulk spin currents driven by an external electrical field for systems 
with either intrinisc or extrinsic spin-orbit interaction. We do not discuss the more 
complicated case when both terms are present. 

When an electrical field is applied and the electrons are accelerated, the pre- 
cession described in Sec. 13.11 was considered. Because the initial spin density for a 
given direction of k is proportional to a, but the non-adiabatic correction is pro- 
portional to 1/a, the spin-orbit coupling constant cancels. Initially, it was believed 
that a small concentration of impurities has no effect and a "universal" spin Hall 
conductivity cr^^ = e/Anh was predicted (Sinova et al, 2004). However, it turns 
out that when the impurities are properly taken into account, the vertex correction 
cancels the bubble term, see Fig. |31 Thus, the dc conductivity vanishes (Inoue 
et al, 2004; Raimondi and Schwab, 2005), 

a^" = 0, (30) 

which was confirmed in numerical calculations (Sheng, Sheng, Weng and Haldane, 
2005). Only when the ac conductivity <t^^ (w) is considered, in the regime l/r ^ 
u) <^b/h the universal value is recovered (Mishchenko et al, 2004). 

That there are no bulk spin Hall currents can be unterstood by the following 
argument due to Dimitrova (2005). Using the Heisenberg equation and for parabolic 
bands, one finds the identity ddy/dt — — i^ma/f?^ {^zVy + ^y^z) for single particle 
operators (Burkov et al, 2004; Erlingsson et al, 2005). For a homogeneous system. 
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Figure 3: Diagram for spin Hall conductivity. The Kubo formula for a is pro- 
portional to Ty {{jyGiij^GA)), with charge current operator and retarded and 
advanced Green's functions Gji/^, and where ((•)) includes averaging over impurity 
configuration. In lowest order in l/k-pi, we can neglect crossed impurity lines and 
p.SH jg given by the diagram shown in (a). Here, the full lines symbolize the renor- 
malized Green's functions including self-energy. In (b), the vertex renormalization 
due to impurity scatterings (connected by dashed line) is defined recursively. When 
only the first term of (b) is taken, we get the bubble contribution to cr^^; when all 
terms are summed, this leads to the additional (ladder) vertex correction. 



one then takes the expecation value of this identity and finds that jy oc dSz/dt. 
When we consider dc properties, we must assume that the system is in a stationary 
state (i.e., we need impurity scattering). Then, the spin polarization Sz is constant 
and thus jy — 0. These arguments, as well as spin current cancellation in a magnetic 
field in absence of scatterers (Rashba, 2004) show that the cancellation is an intrinsic 
property the of free electron Hamiltonian Ha and is not related to any specific 
property of the scatterers. Furthermore, Chalaev and Loss (2005) find that the weak 
localization contribution to jy vanishes, and show more generally that jy vanishes 
even if both Ha and Hf} are present. Grimaldi et al. (2006) find vanishing jy for 
arbitrary values of akp/Ep. 

This cancellation is special to fc-linear spin-orbit interaction (for nonparabolic 
bands, there can be a small finite contribution, proportional to (Krotkov and 
Das Sarma, 2006)). For example, for two-dimensional hole systems, the coupling 
Ha,h is cubic in k [Equation (|lip ]. Then, if isotropic scattering is assumed, the 
vertex correction vanishes, and a "universal" value cr^^ — Se/Anh is found in the 
clean limit br ^ h (Murakami, 2004; Schliemann and Loss, 2005). Quite generally, 
for a 2D system where the spin-orbit field b (k) winds N ^ 1 times around a circle 
in the xy plane when k moves once around the Fermi circle (i.e., A'^ = 3 for Ha,h), 
a universal value 

is found in the clean limit and for isotropic scattering (Shytov et al, 2006; Khaetskii, 
2006). 

Also for the fc^-Dresselhaus couplings, i?D,3d and i?D,2d [Equations ^ and ©J, 
the vertex corrections vanish for isotropic scattering. This leads to a finite spin Hall 
conductivity for two dimensional systems when iJo. 2d is included (Mal'shukov and 
Chao, 2005). Similarly, a^'^ is finite for -ffD,3d (Bernevig and Zhang, 2004). 
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The results cited above are only valid for isotropic scattering, except Equation 
([SO]), which holds more generally. Recently, more general descriptions using kinetic 
equation (cf. Sec. 13. 4p allowed to include arbitrary angular dependence of impurity 
scattering (Shytov et ai, 2006; Khaetskii, 2006). It turns out that cr^^ significantly 
depends on the shape of the scattering potential and does not reduce to a simple 
form in general. For example, in the clean limit br h and in the regime of small 
angle scattering (but still for a typical momentum transfer q > b/vp, i.e., not too 
small angles), one finds 



S(^)(^-^-) 



where C describes the non-parobalicity of the band, oc k^'^'', and for |b| c>c 
(Shytov et al, 2006). For example, taking C = and Ha^h, with iV = = 3 , we see 
that the sign of in Equation [Equation ([32| ] is opposite to the case of isotropic 
scattering [Equation (|3ip ]. Similarly, Liu and Lei (2005) found in a numerical study 
of system with spin-orbit coupling Ha,h that cr^^ strongly depends on the type of 
the scattering potential. Therefore, the spin Hall conductivity is not a universal 
quantity, as its numerical prefactor and its sign depend on sample parameters. On 
the other hand, for clean systems with N ^ 1, the order of magnitude is consistently 
|aS"| - e/A-Kh. 

The above results are valid for weak spin orbit coupling, ka <C k-p. Conversely, 
there are also materials with strong spin-orbit coupling, as it was recently found 
in Bi/Ag (111) and Pb/Ag (111) surface alloys (Ast et al, 2005). This motivated 
Grimaldi et al. (2006) to generalize the theory; however, relying on an extensive 
numerical procedure. 

The extrinsic contribution Hcxt.c for electrons also leads to spin currents 
(D'yakonov and Perel', 1971; Hirsch, 1999). These currents are often evaluated only 
for isotropic impurity scattering (Zhang, 2000; Shchclushkin and Brataas, 2005). 
Assuming absence of intrinisc spin-orbit interaction, for aribitrary angular 
dependence of scattering, the extrinsic spin Hall conductivity equals (Engel 
et al, 2005) 

(T^H - -^a,, + 2nA|, (33) 

where the first term is due to skew scattering (see Sec. l3.2p . The second term due to 
the side-jump mechanism (Sec. l3.3p : as this mechanism goes beyond transport equa- 
tion, this term has to be evaluated separately. Here, a^x is the electrical conductivity 
and we defined the transport skewness 

^ Jdni{^) [1- cos ^) ^ ' 

that depends on the structure of the scattering center and on the Fermi energy, 
and /, S are defined below in Equation pT|) . For screened Coulomb scatterers. 
Equation ((33|) can be evaluated without any free parametrers (Engel et al, 2005) and 
the resulting absolute value of spin current is in quantitative agreement (within error 
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bars) with the observation by Kato et al. (2004b) in GaAs and seems comparable with 
the data by Stern et al. (2006) in ZnSe — implying that the observed spin currents are 
due to the extrinsic effect. Note that in Equation ([55)) the skew scattering and the 
side jump contributions have opposite signs. The skew scattering term dominates in 
standard transport theory where one expands in h/ Et where E is, a, typical electron 
energy; however, for dirty samples both terms can be of comparable magnitude. 
Stern et al. (2006) found that the measured cr^^ in ZnSe has the sign of the skew 
scattering contribution and that the same is likely to be the case for the cr^^ observed 
by Kato et al. (2004b). Finally, assuming short-range scatterers, Tse and Das Sarma 
(2006a) found the same order of magnitude for extrinsic spin currents. However; 
later they concluded that intrinsic spin-orbit coupling can cancel skew-scattering 
and reduce side-jump contributions to cr^^ (Tse and Das Sarma, 2006b). 

Remarkably, despite the fact that the side jump term in Equation (j33p was de- 
rived by including electron dynamics during the scattering event, it does not contain 
any factors related to the scattering probability — it only depends on the coupling 
constant A, which is an intrinsic property of the material and is directly related to 
a Berry connection through the spin-orbit contribution to the operator of electron 
coordinate, cf. Sec. 12.61 Thus, while in this review and commonly in the literature 
the side-jump contribution is considered as extrinsic, it is clear that the distinction 
between intrinsic and extrinsic is somewhat arbitrary for this contribution. 

4.4 Anomalous Hall effect and its relation to spin Hall effect 

In the anomalous Hall effect ( AHE) , equilibrium polarization of a ferromagnet com- 
bined with spin-orbit interaction leads to electrical Hall currents transverse to an 
applied field. The theory of AHE has a long history and reveals many problems 
typical of spin transport in media with spin-orbit coupling, including the competing 
mechanisms of spin-orbit scattering by impurities and the role of intrinsic spin pre- 
cession; for reviews see (Nozieres and Lewiner, 1973; Crepieux and Bruno, 2001; Na- 
gaosa, 2006) . For non-interacting electrons and negligible spin relaxation, AHE and 
SHE are closely related; this is true for extrinsic spin-orbit interaction because A is 
small and spin relaxation is of order A'^ (Elliott, 1954; Yafet, 1963). In the SHE, 
we can decompose the spin currents j as a difference in particle currents of two 
spin species with polarizations ±e^. Regarding these species separately, each carries 
the anomalous Hall current J^'^ of a system with spins fully aligned along the ±£^ 
direction and with density txah = \n, because we consider the SHE in non-magnetic 
media, where electrons are unpolarized in equilibrium. We can express the spin Hall 
current as (Engel et al, 2005) 



This relation allows to make use of the extensive literature on the AHE to gain 
further insights into mechanisms of the SHE, at least on its extrinsic part. 




(35) 



21 



4.5 Spin accumulation and transport at boundaries 

For only extrinsic spin-orbit coupling, because the spin relaxation is negligible for 
small A, the spin is almost a conserved quantity and thus spin density and spin 
current satisfy a continuity equation with a small relaxation term. In this case, bulk 
spin currents will produce a spin polarization at the edge of the sample, i.e., spin 
currents and spin accumulation are directly related (D'yakonov and Perel', 1971; 
Hirsch, 1999; Zhang, 2000). The polarization at a y = edge is Sz = {h/2)sz = 
{h/2) ^jTs/Dsjy, with spin relaxation time Tg and spin-diffusion coefficient Ds (which 
is identical to the electron diffusion coefRcien D in the absence of electron-elecon 
interaction) . 

For intrinsic spin-orbit interaction, it is not clear whether any general relation 
exists between spin accumulation and bulk spin currents, but spin accumulation can 
be studied directly. We discuss here the situation of a semi-infinite two-dimensional 
electron system, with Rashba coupling, located in the upper half-plane {y > 0). 
We assume a uniform applied electric field parallel to the x-axis, and we consider 
the diffusive limit, where the spin diffusion length is large compared to the mean 
free path, so Eqs. ((25)) and (|26|) apply far from the boundary. The spin density 
near the edge will depend on the boundary conditions to the diffusion equations at 
y = 0, and these will depend, in turn, on the boundary conditions of the microscopic 
Hamiltonian, as has been discussed in the various articles cited in the last paragraph 
of Sec. E31 

In the case of an ideal reflecting boundary, the spin density s(jj) is found to be 
constant, and the same as in the bulk, right up to the edge (Bleibaum, 2006). Thus, 
one finds Sz = 0, while Sy has a value proportional to the charge current and to the 
Rashba coupling constant a. By contrast, if there is strong spin-orbit scattering at 
the boundary, all components of s should vanish there. In this case, the coupled 
diffusion equations predict that for y > 0, there will be non-zero values of both Sz 
and Sy, with oscillating behavior, in a region near the edge whose width is of the 
order of spin-precession length k^^ = fi^ /ma, which is about the Dyakonov-Perel 
spin-diffusion length y/DgTs (Rashba, 2006). 

What happens if the boundary at y = is partially or completely transmitting, 
and there is a second conductor in the region y < which has no spin orbit coupling? 
As noted by Adagideh and Bauer (2005), one should expect, in general, to find non- 
zero oscillatory values of both Sz and Sy in the Rashba conductor near the boundary, 
and injection of spin into the non-spin-orbit material. However, it was found by 
Tserkovnyak et al. (2006) that this actually will not happen in the simplest case 
that one might consider: a boundary between a pair of two-dimensional systems 
with identical properties except for different values of a. There will quite generally 
be a discontinuity in the spin densities at a lateral boundary between systems with 
different Rashba coupling. If the electron mobility is a constant across the boundary, 
the discontinuity turns out to be equal to the difference between the bulk spin- 
densities of the systems. Then, on each side of the boundary one finds Sz — 0, while 
Sy is the same as the respective bulk value. There will thus be no spin injection if 
the second two-dimensional electron system has a = 0. 

Mal'shukov et al. (2005) have argued that there should be spin accumulation 



22 



near a reflecting edge, in the diffusive case, when the cubic Dresselhaus coupling 
Ht), 2d is important. In the opposite hmit of balhstic transport and near a sharp 
specular edge, Usaj and Balseiro (2005) have found spin magnetization due to Ha 
that oscillates rapidly with a period of fcp ^ and shows beating on a length scale of 

Ha ■ 

There are also numerical approaches analyzing the edge spin accumulation. No- 
mura et al. (2005) simulated a two-dimensional hole system, using the coupling 
Ha,h [Equation (jlip j. They found a spin accumulation that was consistent with the 
experiments by Wundcrlich et al. (2005). 

4.6 Mesoscopic systems and spin interferometers 

So far we have not considered interference effects — e.g., in Sec. 13.41 we used an 
expansion in lowest order 1 / fcpf which does not include interference between electron 
propagation paths that follow different trajectories. To include such coherent effects 
in systems with impurities, one needs to consider the next order in l/kpi- the 
weak localization corrections. Alternatively, one can consider clean systems that 
are ballistic on length scales of the device. Spin interference effects in mesoscopic 
systems open up a new set of technical possibilities, e.g., in rings or ring- like arrays 
with spin-orbit interaction one can use Berry phase and Aharonov-Casher phase 
effects to study a variety of phenomena. In the presence of an applied magnetic 
field, spin-orbit effects can modify the Aharonov-Bohm or Altshuler-Aronov-Spivak 
oscillations in the electrical conductance. For theoretical discussions, see (Aronov 
and Lyanda-Geller, 1993; Bulgakov et al, 1999; Engel and Loss, 2000; Frustaglia 
et al, 2001; Koga et al, 2004; Aeberhard et al, 2005). For experimental results, see 
(Morpurgo et al, 1998; Yau et al, 2002; Yang et al, 2004; Koga et al, 2006; Koenig 
et al, 2006; Bagraev et al, 2006). 

For practical applications, it is not only important to generate nonequilibrium 
spin polarization in media with spin orbit coupling, but also inject spin currents 
produced by such populations into "normal" conductors, i.e., conductors with negli- 
gible spin orbit coupling. In normal conductors spin is conserved and spin currents 
are well defined. Spin injection can be achieved using spin-orbit coupling, even in 
devices without magnetic fields and without ferromagnetic components. Propos- 
als for such devices, in the mesoscopic regime, have been made by (Kiselev and 
Kim, 2003; Shekhter et al, 2005; Souma and Nikohc, 2005; Eto et al, 2005; Silve- 
strov and Mishchenko, 2006). 

Generally, spin interference devices make use of the intrinsic and extrinsic spin- 
orbit couplings presented in Sec. [2] and the spin transport mechanisms discussed 
in Sec. [3] are important. However, we do not present more concrete descriptions of 
interference effects or details of microscopic structures; this could be done by solving 
the Schrodinger equation analytically, by simulating it numerically, or by using weak 
localization calculations. On the other hand, we can assess the length scales on 
which spin precession effects are expected: both diffusion equation (Sec. 13. 5p and 
its solution near boundaries (Sec. 14. 5|) indicate spin precession length £a = ^/ka 
as a characteristic length for spin distributions. Furthermore, for clean systems, 
responses to an inhomogeneous field diverge at the wave vector q = 2ka of the field 
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(Rashba, 2005). This "breakdown" suggests the length as an optimal size for 
achieving large spin polarizations. For a more generic Hamiltonian Hint, this scale 
can be estimated as icS ~ fi^kF/iTilbl, establishing a "mesoscopic" scale at which 
one can expect largest static spin responses to electric fields. Because ieS is also of 
the order of the Dyakonov-Perel spin diffusion length, this estimate seems applicable 
both to the ballistic and diffusive regimes. 

5 Spin Hall effect due to edge states in insulators 

In the previous parts of Sec. [U we discussed spin currents jy driven by electric field 
Ex, their relevance to spin transport and spin accumulation, and also the techniques 
for calculating conductivities [Equation (|29p ]. Even when these nondiagonal 
components of the tensor jj were not directly influenced by dissipation, they were 
calculated for ordinary metallic conductors whose longitudinal electric conductivity 
a^x was controlled by electron scattering, hence, electron transport in the bulk was 
dissipative. More recently, Murakami et al. (2004) proposed that some centrosym- 
metric 3D systems possess properties of "spin insulators." These are media with 
gapped electron spectra and zero bulk electrical conductivities Uxx but finite and 
dissipationless spin conductivities ct^^. 

The basic idea is as follows. A set of electron bands that in absence of SO coupling 
belongs to orbital momentum L, in presence of SO coupling is described by the total 
angular momentum J = L+S. When some of the bands belonging to the J multiplet 
are filled, while different bands of the same multiplet are empty and are separated 
from the filled bands by a gap, all filled bands contribute to spin current. Uniaxially 
strained zero-gap semiconductors a-Sn and HgTe, and narrow-gap semiconductors 
of PbTe type were proposed as model systems. Spin conductivity cr^^ is large in 
these materials: it is about e/ha in 3D, a being the lattice constant. In 2D, cr^^ 
is quantized when the Fermi level is inside the gap (Qi et al, 2006; Onoda and 
Nagaosa, 2005). However, for reasons similar to those discussed in Sec. 14.11 the 
relation of this cr^^ to spin transport is not obvious and was already questioned 
(Kane and Mele, 2005a; Kane and Mele, 2005b). 

A different concept of spin transport in an insulating phase has been developed 
by Kane and Mele as applied to graphene (Kane and Mele, 2005a; Kane and Mele, 
2005b). It is based on the Haldane model of quantum Hall effect (QHE) with 
spinless fermions under the conditions of zero total magnetic flux across the unit 
cell (Haldane, 1988). It has been emphasized (Onoda and Nagaosa, 2005; Kane and 
Mele, 2005a) that this model differers fundamentally from the model by Murakami 
et al. (2004), in particular in the properties of edge channels. Their crucial role for 
the QHE has been clarified by Halperin (1982), and they play a similar role in the 
physics of spin Hall effect in graphene. In what follows, we consider properties of 
graphene in more detail. The graphene model is not only of conceptual interest but 
is also attractive because of the very recent experimental achievements in measuring 
electron transport in graphene (Novoselov et al, 2004; Zhang et al, 2005). 

Graphene is a monoatomic layer of graphite. Its honeycomb 2D lattice is shown 
in Fig. m The elementary cell includes two atoms shown as A and B. The phase 
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Figure 4: Schematics of honeycomb lattice of graphene. The hexagon in the center 
is an elementary cell containing two carbon atoms that belong to two sublattices. 
Atoms of these sublattices are marked as A and B and are shown by empty and filled 
circles, respectively. Brillouin zone is shown by a dashed line. K and K' indicate 
nonequivalent corners of the zone where the gap opens. 



diagram of graphene can be understood from the tight-binding Hamiltonian (Kane 
and Mele, 2005b) 

^ = i^4<^3 + -jn^so X! ^yc+cr^cj + i\R ^ c+(cr X dij)^Cj + A„ ^ C^c+Q , 

(36) 

where Ci are annihilation operators at lattice cites i, spin indeces in them being 
suppressed. The first term is the nearest neighbor hopping term between two 
sublattices. For the following, of the critical importance is the second term with 
Vij — (2/\/3)(di X d2)z = ±1 that describes second neighbor hopping. Here di and 
d2 are unit vectors along two bonds that an electron traverses when going from the 
site i to the site i. The cross product of di and d2 produces a screw that in the Hal- 
dane model of spinless fermions couples them to inhomogeneous magnetic flux, while 
in the present model it couples the orbital motion of an electron to Pauli matrix ■ 
Hence, Xso is a coupling constant of a mirror symmetric, z —z, spin-orbit inter- 
action. The third term is a nearest neighbor Rashba term, d^j being a unit vector 
in the direction connecting i and j nodes. It explicitly violates z — > ~z symmetry 
and originates from the coupling to the substrate or from an external electric field. 
The fourth term is a staggered sublattice potential with taking values = ±1 for 
A and B lattice sites. It vanishes for graphene but would be present for a similar 
boron nitride BN film. Including this term is a clue for explaining the difference 
between the "quantum spin Hall" (QSH) phase and a usual insulator (Kane and 
Mele, 2005a; Kane and Mele, 2005b). 

The remarkable spin properties of graphene are seen from the one-dimensional 
projection of its energy spectrum. Fig. [5l found by solving the Hamiltonian of 
Eq. (|36p in the geometry of a strip with finite extent in the y-direction (defined 
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Figure 5: Energy bands of a one-dimensional strip of graphene with "zig-zag" edges, 
i.e., the strip has finite extent along the y-direction shown in Fig. [H Narrow gaps 
in the bulk spectrum are achieved at the K- and K'-vertices of the Brillouin zone. 
Branches of the spectrum originating from the edges of bulk continua show energies 
of edge states. Edge states at a given edge of the strip cross at ka = n, a being 
the lattice spacing, (a) QSH phase for A„ = O.lt; edge states at a given edge of 
the strip cross at ka = n, a being the lattice spacing, (b) A normal insulating 
phase for = OAt. In both cases Xso = 0.06t and Xr ~ O.Q5t. The inset shows 
the phase diagram in the Xy-Xn-plane for < Xso ^ t. Reprinted figure with 
permission from [C.L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 146802 (2005)]. 
[Copyright (2005) by the American Physical Society., 

in Fig. H]), i.e., having "zig-zag" edges aligned along the x-direction. The spec- 
trum comprises four energy bands, of which the two lower bands are occupied; each 
bulk band is two-fold spin degenerate. Narrow gaps at K and K' open because 
of Aso,A„ ^ — for Xso,^v = 0, electrons possess a k-lmeax spectrum of Dirac 
fermions, e{k) = hc*k. In addition to bulk states, there are edge states connecting 
K and K' bulk continua. Their topology in the panels (a) and (b) is rather different. 

In the panel (a) drawn for a small X^, edge states traverse the gap. For each edge 
of the strip, there are two such states. They are Kramers conjugate and propagate 
in opposite directions. This behavior reflects unusual cross-symmetry of bulk states 
that manifests itself in the opposite signs of the gap function at K and K' points. 
The small- A^ phase has been dubbed as QSH-phase by Kane and Mele (2005a). It 
is the distinctive property of this phase that at any energy inside the gap there is 
one pair of edge modes (more generally, an odd number of such pairs). 

When Xr — 0, cr^ is conserved, and the pattern of dissipationless spin transport 
become especially simple. Each of independent subsystems of cr^ =t and cr^ —I elec- 
trons is equivalent to Haldane spinless fermions (Haldane, 1988). One pair of such 
"spin filtered" states propagates along each edge. The states with opposite az po- 
larizations propagate in opposite directions. Because they form a Kramers doublet, 
potential backscattering is forbidden and transport is dissipationless. This model 
predicts two-terminal electric conductivity G = 2e^/h. Propagation of charge cur- 
rent through edge states results in antisymmetric spin accumulation at these edges. 
In four-terminal geometry, spin currents flow between adjacent contacts, and related 
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spin conductances are quantized; when normalized on a number of transported spins, 
= ±e/27rS. For Xr ^ 0, however, ct^ is not conserved. Nevertheless, spin currents 
persist (if Xr is small, see below) but they are no longer exactly quantized. An early 
argument that the spin Hall conductance can be quantized was given by Froehlich 
et al. (1995) when considering incompressible 2D systems. 

Panel (b) of Fig. O drawn for a larger staggered potential A„, shows properties 
of a normal narrow gap insulator. The gap function has the same sign at K and K', 
as a result, one pair of edge states runs between two conduction band valleys and 
the second pair between two valence band valleys. For some boundary conditions 
at the strip edges, edge states can penetrate the gap. However, there is always an 
even number of Kramers conjugate pairs of edge states at any given energy inside 
the gap, hence, backscattering is no longer forbidden. Therefore, it is the topology 
of edge states that defines the difference between the QSH and insulating phases in 
a simple and concise form. 

The QSH phase is formed due to the bulk spin-orbit coupling Xso- Increasing 
asymmetric Xr or staggered A^, potentials destroy it when they become large enough. 
A phase diagram of the competing phases, QSH phase and a normal insulator, 
is shown in the inset to Fig. [S] The QSH phase exists inside an ovaloid in the 
Xy/ Xso~^rI ^SO plane. Outside it, graphene shows properties of a normal narrow 
gap insulator. 

Another factor suppressing the gap and spin conductivity is electron scattering in 
the bulk. Its effect has been investigated numerically by Sheng et al. (Sheng, Sheng, 
Ting and Haldane, 2005) for a four-probe spin Hall setup by using the Landauer- 
Biittiker formula (Landauer, 1988; Biittiker, 1988); their spin conductivity tr^^ de- 
scribes real spin transport. Disorder was modeled as eicfci with ei randomly 
distributed in the interval [—W/2, W/2]. They found that remains within a few 
percent of its quantized value when W < t and the Fermi level stays inside the gap; 
and cr^^ drops fast with increasing W for W ^ IM. Under the same conditions, 
(7^^ remains stable for Xr < 0.2t. These results allow to establish the parameter 
range inside which edge spin channels remain robust and carry dissipationless and 
nearly quantized spin currents. Inside this range, there exist close analogy between 
QSH and QHE systems. 

Currently, there is no direct experimental indications of the spin gap in graphene. 
A crude theoretical estimate of it by Kane and Mele (2005a) results in the gap 
2Aso ~ 2.4 K what is in a reasonable agreement with different data (Brandt 
et al, 1988). However, more recent calculations indicate that the actual value of 
Aso is considerably smaller. Meantime, estimates of H/t based on transport data 
result in H/t > 25 K for typical mobilities of « 10,000 cm^ s^-^, t being the 
momentum relaxation time. Comparison of these data can easily explain suppression 
of spin-polarized transport through edge channels by disorder in the samples that 
are currently available. An independent mechanism of suppression is the Xr con- 
stant that develops when electron concentration is controlled by a gate, and the ratio 
Xr/Xso is unknown. Unfortunately, all estimates are crude because electron trans- 
port in graphene is still not understood. Novoselov et al. (2005) recently reported the 
minimum metallic conductivity 4e^/ft. when the Fermi level passes through the conic 
point of the spectrum; it is nearly independent of the mobility /i. A nonuniversal 
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and even larger conductivity of about 6e^//i was reported by (Zhang et al, 2005). 
The closest theoretical value 2e^/h comes from the spin channel model (Kane and 
Mele, 2005a), followed by 4e^/7r/i and 2e^ /nh found from different models of the 
bulk transport of Dirac fermions (Zicglcr, 1998; Shon and Ando, 1998). Because 
there exists a region of parameter values where spin transport through edge chan- 
nels is robust, honeycomb lattices lithographically produced from semiconductors 
with strong spin-orbit coupling may also be of interest (Zheng and Ando, 2002). 

Abanin et al. (2006) argued that in the QHE regime, the exchange-enhanced 
gap for chiral edge modes, originating from Zeeman splitting, may be as large as 
100 K. Another system where edge spin channels may play a role was proposed by 
Bernevig and Zhang (2006); it includes parabolic confinement in conjunction with 
inhomogeneous shear deformation. More recently, Fu and Kane (2006) proposed 
Bii-ajSba; semiconductor alloys and a-Sn and HgTc imdcr uniaxial strain as materials 
that satisfy the necessary symmetry requirements for the QSH phase and that are 
expected to have a large Aso- Independently, Bernevig et al. (2006) showed that 
HgTe/CdTe quantum wells are also good candidates for the QSH phase. 
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